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Tbi*  Memorandum  daseribaa  two  computer  programa  which  bar*  baan  develop  ad  aa 
data  preparation  toola  for  tha  three-dimensional  boundary  integral  aquation 
(BIE3D)  conputar  program.  Tha  firat  ia  a  three-dimensional  mesh  generation 
conputar  program  which  ia  baaed  on  a  semi-automatic  subdivision  technique  similar 
to  a  two— dimensional  finite  elanent  aash  generation  procedure.  The  progran  ia 
written  in  standard  FORTRAN  and  it  ia  now  operational  in  the  Materials  and 
Structures  Departnantal  TAX  computer.  The  second  la  a  conputar  progran  for  plot¬ 
ting  the  generated  mesh.  Tha  plotting  program  is  used  to  drew  a  three-diaanaional 
wiew  of  the  mesh  of  a  problem.  In  particular,  a  generated  neah  can  be  drawn  so 
that  its  topology  can  be  checked  visually  for  anomalous  points  and  incompatibility 
of  elements.  The  program  is  written  in  BASIC  and  it  runs  interactively  on  a 
Hewlett-Packard  HP26477  graphics  terminal . 

Both  computer  programs  are  described  with  examples.  Instructions  and 
criptions  for  both  input  and  output  data  file  are  presented  in  each  case. 
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i  urrraouciHM 

TUi  Memorandai  describes  two  wpoMt  pcogfs  tAiek  km  boon  developed  *<  data 
preparation  tool*  for  the  three-dimensional  boundary  integral  equation  (BH3D)  computer 
program1 .  These  progrm*  will  help  IB3D  nears  to  perform  their  teaks  wore  efficiently 
and  effectively.  In  particular,  tiae  savins  is  the  biggest  advantage  of  using  these  pro- 
grans.  Tbs  programs  are  called  MMDD  and  SIAR3D. 

MESH3D  ia  a  three-dimensional  assh  generation  computer  program  for  use  with  the 
three-dimensional  boundary  integral  aquation*  computer  program  OZX30) ’ .  It  is  written 
ia  standard  FOKTUI  and  it  is  available  ia  the  Materials  and  1  true  turns  Departmental  Viz 
computer. 

MAM3D  is  a  computer  program  that  can  be  used  to  plot  a  generated  three-dimensional 
mash  on  a  graphics  visual  display  or  on  paper.  The  program  is  written  in  BASIC  and  it 
can  be  operated  interactively  by  the  users  on  a  Bewlett-Facksrd  IP2A47T  graphics  tsrninal . 

These  prograns  are  now  available  to  all  BZZ3D  users  as  an  option.  The  users  can 
either  create  a  HZ3D  input  data  fila  ia  tha  conventional  way,  as  described  fully  ia  a 
previous  Memorandum1,  or  use  the  IWIH)  program  as  described  here  ia  section  2.  A  graphic 
output  of  the  generated  mesh  can  be  Obtained  using  the  1MB  program,  as  described  in 
section  3. 

The  main  aim  of  this  Wsmorandiai  is  to  lay  dona  tha  guide  lines  of  the  operating 
procedures  of  tbs  prograns.  Step-by-step  Instructions  are  given  to  show  how  an  input  and 
output  data  files  can  be  prepared  for  both  of  tbe  programs.  Furthermore,  test  cases  for 
both  of  ths  ytcgtms  era  illustrated  with  smanples  and  descriptions. 

Za  order  to  redoes  the  amsuat  of  time  involved  in  data  preparation  for  a  boundary 
integral  aquation  (BIB)  analysis,  rsmputsr  aided  design  (CAD)  software  soch  as  mesh 
generation  computer  programs  can  ha  employed.  The  use  of  soch  programs  not  only  saves 
tbs  but  also  eliminates  the  haaa  errors  which  occur  during  manual  preparation  of  the 
nodal  co-ordinates  and  element  topologies.  Mesh  generation  prograss  hsve  been  well  docu¬ 
mented  in  the  finite  element  nathsd  (FBO.  These  prograns  are  generally  of  two  types2’*: 

(a)  A  digitising  tablet  coupled  with  an  interactive  graphics  system  may  be  alloyed  to 
define  and  produce  tbs  geometric  data.  A  light  pen  is  usually  used  to  digitise  tbe 
structural  outlines  of  a  preliminary  assh.  Tbs  program  nay  than  allow  for  the  choice  of 
element  types  and  ths  fineness  eC  subdivision  of  elements.  Vith  this  information,  tbe 
mesh  generation  program  will  generate  a  mash  antemsticelly  and  output  it  onto  a  graphics 
visual  display  or  a  plotter.  However  assh  generation  algeridms  of  this  type  are  highly 
dependent  on  the  neks  and  medal  of  the  digitiser  and  normally  require  specialised  system 
software. 

(I)  A  semi  an rums tic  approach  may  be  used,  where  tha  structural  outline  is  based  on  a 
faw  large  blocks  which  are  defined  as  input  data.  Blsnsets  are  than  generated 


*  This  method  is  also  referred  to  as  the  boundary  llsaont  Method  (BBM>. 
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— t— tlcally  by  a  ■UlTiiioo  pneMi.  Xa  this  subdivision  process,  tha  generated 
■lant  alaaa  aaai  to  la  qaeUlad  aa  a  aaparata  ut  of  Input  data.  Tbia  approach  ia 
described  fa  da tall  fa  tha  following  sections  and  It  la  tha  aathod  uaad  fa  tha  praaaat 
study.  d  fa  dtsaaataaal  ffaita  alaaaat  aaah  gaaaratloa  routine*  which  ia  baaad  on  thia 
approach  haa  haaa  adoptal  aal  aodtflad  for  generating  a  aaah  on  a  two-diaaaaional  surface, 
aa  required  far  a  throw  fltaanalnnal  HX  calculation. 

2.1  Tha  aaah-aaaof  tfon  esunfrr  proaraa  0C8H3B) 

Thia  rhraa  diaaaaiaoal  aaah  ganararlrm  program  haa  hoan  developed  jointly  with  tha 
Ooparflaaat  of  tagisaaring  Materials,  University  of  goutha^tou  and  ia  now  inatallad  in 
tha  Matoriala  and  Structures  hoparfantal  VAX  cowputor.  Tha  protraa  CKXSH3D)  ia  baaad  on 
autoaatfc  alaaaat  subdivision  of  a  faw  largo  atructural  hlocka  which  ara  defined  aa  input 
data!  each  block  la  auHtridad  fate  alawaata  according  to  weighting  factora  provided  aa 
input  data. 

«*— Q"  will  calculata  all  tha  nodal  co-ordinate*  for  all  tha  auhdividod  alaaanta. 
The  nodal  point  adhering  eonvaation  Celockwiaa  or  oaticloclarlaa) ,  geaoratad  for  tha  aaah, 
will  he  tha  sons  aa  tha  input  data  iwhartng  eoavention  for  tha  atructural  blocks. 

2.2  gawaerlc  ranraaantation  of  a  atructural  block 

In  tha  netted  of  Big,  tte  boundary  contour  ia  the  aurfaca  of  a  three-dinaaaioaal 
body  and  ia  alauya  rapraaantad  By  aurfaca  alauanta  taxaad  boundary  eluants.  Tha  type  of 
boundary  eluent  uaad  ia  RX3D  ia  tha  eight-node  loop  arena  trie  quadrilateral  alanent,  aa 
ehoun  in  Tig  1 .  MESH3S  uaao  tte  aight-noda  quadrilateral  to  represent  a  atructural  block 
in  tha  initial  ataga  of  tte  aaah  generation  process.  The  nathanatical  representation  of 
a  structural  block  to  based  on  tha  global  cartesian  co-ordinate  eye  tan  where  each  of  the 
nodes  *B  to  appraadnatad  by  interpolation  forwula  which  involves  finite  olswaat  shape 
function#  ao  shown  balowt  . 

*<£,«o  -  Y  0) 

whore  ja  arc  tha  nodal  co-ordinataa  and  M°CC,n)  arc  the  chape  functions*  for  a  two- 
dineaaional  finite  elenent  analysis  x  »  Cx.y,*)  . 

This  intarpolation  io  equivalent  to  eh#  eranafornaeion  of  a  boundary  elenent  to  a 
piano  oquare  in  a  local  co-ordinata  ayatan  (C.n),  aa  shown  in  Fig  1.  For  axanpla,  in 
Vig  I,  the  lino  X  •  ♦!  will  wap  into  tha  global  co-ordinaea  eye ten  such  that  nodes  5. 
4,  3,  will  lis  on  the  sons  curve i  this  ia  aa  appralastioa  to  the  actual  ahapa  of  the 
edge  of  tha  boundary  alanent.  Other  lines  C.n  •  constant  will  sinilerly  nap  into 
curves  ia  the  global  eo-ord lasts  eye  ten  to  approninaea  the  boundary  alanent,  which  is 
being  paranetrically  rapraaantad  by  f,n  . 

2.3  Tha  subdivision  process 

Any  one  aight-noda  quadrilateral  alanent  nay  be  subdivided  into  a  outer  of  ualler 
alanent a  and  tha  else a  of  the  subdivided  alanent*  are  prede t trained  by  waigbeing  factors 
specified  aa  input  data. 
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When  an  element  la  subdlvi4ed  astra  nodal  will  ba  craatad  in  order  Co  ratcora  tha 
coupatibility  of  tha  eight-node  alasasca.  Tor  esople,  in  subdividing  a  tbraa-dinaniional 
tody,  such  aa  a  onto  in  Fig  2,  tto  anbdiriiion  of  ona  of  tto  aliments  on  ona  aida  of  tha 
caba  casaaa  incompatibility  of  nodal  poises  with  tto  elements  adjaeant  to  it.  Further 
additional  nodal  point a,  in  all  a  total  of  twelve,  suit  ba  introduced  in  ordar  to  rastora 
tto  compatibility  of  tto  elements.  In  fact,  tto  subdivision  procaaa  not  only  subdivide* 
a  largo  block  into  mmller  alaaonta  tot  alao  ensures  all  tto  elements  and  all  tto  nodal 
points  ara  conpatibla  with  sack  other;  rapaatad  canon  nodal  at  boundary  intarfacaa  art 
altoiaatad  Caaa  aaetlon  2.4).  Sn  ordar  to  achieve  this,  tto  aubdlTiaion  procaaa  requires 
a  act  of  input  data  that  itself  ia  consistent  and  conpatibla.  For  instance,  for  each 
structural  block  it  is  necessary  to  specify  tto  following: 

(a)  Tto  ""»*»»  of  subdivisions  to  ba  node  in  both  C  and  n  directions  ara  stored  in 
KDITX  and  asZTT  respectively.  For  aacaaapla,  if  DIB  “  3  and  IDXVT  ■  2  then  tto  block 
will  to  subdivided  into  three  divisions  ia  tto  (  direction  and  two  divisions  ia  tto  n 
direction  respectively,  da  a  total  of  six  elaaents  will  be  generated  as  shown  in  Fig  3. 
Tto  dir  action  of  tto  sobdiviaiona  nut  to  clearly  defined  where  there  ara  too  or  aore 
blocks  linked  together  because  tto  compatibility  at  tto  bloek-to-block  intarfacaa  auat 
be  preserved,  <i*  tto  natar  of  alaaanta  at  tto  adjacent  blocks  anst  be  tto  Sana,  as  indi¬ 
cated  in  Fig  4. 

(b)  weighting  factors  ara  specified  as  input  data,  for  each  individual  block,  to  define 
the  relative  sixes  of  tto  olaanta  within  a  Block.  Tto  weighting  factors  are  denoted  by 
CU{)t  where  X  is  tto  direction  of  tha  weighting  factors.  Ibis  information  is  stored  in 
the  following  formats 

0KHX0om,cn>iv*-J,«Dm)  . 

Similarly,  tto  weighting  factors  ia  the  n  direction  are  stored  ass 

(strrr(ium) ,  (n>m-i  ,mvm  . 

It  should  ba  noted  that  ttasa  weighting  factors  are  used  to  define  tto  relative  sites  of 
tto  elenents  ia  a  particular  easy  for  a  particular  blocks  therefore  their  absolute  values 
do  not  affect  tto  aash  generation  process.  For  aaanpls,  in  order  to  subdivide  a  largo 
structural  block  into  three  elements  tto  weighting  factors  00 ^  nay  be  specified  as: 


10  10  20 


I  1  2 


■o never.  it  is  squally  important  to  specify  all  the  weighting  factors  consistently  for 
all  blocks  adjaeant  to  each  ottor  where  compatibility  of  nodes  at  block-to-bloek  inter¬ 
face  asst  be  reserved.  A  pictorial  illustration  of  compatibility  of  nodes  at  block-to- 
blsck  interfaces  ia  shewn  in  Fig  4, 


lech  local  eo-ordiaatss  C  cal  n  or*  initialised  to  the  value  of  -1 , 
lncronontad  according  to 


«1  -  «i 


0L>1 

"l  "  ”1-1  *  “7“  ' 


i-SDITX  i-HDJVT 

£  <Vt  -  ^  -  £  <Vi 


scion  is  node  over  the 


and  wbors  subscript  1  daaotoa  tbs  ganaral  Inrrnsonf  and  tbs  msaation  is  Bade  orsr  tbs 
total  mber  of  subdivided  elsnenta  la  tbs  specified  direction.  The  real  space 
co-ordinates  are  then  calculated  using  (I). 

2.4  Connectivity  of  individual  blocks 


The  subdivision  process,  as  described  in  the  previous  section,  is  performed  indi¬ 
vidually  for  each  structural  block.  Tbs  generated  alanent  nodal  points  within  each  block 
are  sobered  individually.  Boons ar  alanent  nodal  points  betosen  block-to-block  inter¬ 
faces  will  have  rnfiii  co-ordinate  values  so  they  need  to  be  uniquely  sobered.  This  is 
dona  autoantically  by  se arching  for  all  tbs  points  having  cannon  co-ordinate  velues  and 
replacing  the  nodal  sabers  by  a  single  nodal  noaber.  It  should  be  noted  that  this  pro¬ 
cess  of  allninatiou  of  repeated  nodes  will  in  general  produce  a  new  nodal  sobering  system 
for  tbs  generated  assb,  even  for  the  ease  where  all  the  structural  blocks  ranain  undivided; 
an  asaaple  of  this  is  illustrated  in  Tig  Sabb. 

i  typical  asSh-generation  routine  nay  be  illustrated  using  tbs  siapls  flow  diagran 
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2.3  Input  data 

Subroutine  XH8L0C  reads  In  tha  data  necessary  to  define  a  preliminary  math  which 
may  contain  several  structural  Blocks.  In  the  MESH3D  input  date  file,  apart  from  the 
problaa  title,  all  data  may  be  entered  as  free-fornat. 

Descriptions  of  a  HESB3D  input  data  file  is  given  below  and  in  order  to  aee  bow 
this  works  a  test  case  is  Illustrated  in  section  2.6  with  tha  entire  data  file  of  the 
problem  listed  with  descriptions.  The  users  may  use  this  example  as  a  guide  line  to 
prepare  a  data  file  for  tha  XESH3S  computer  progra. 

For  clarity,  the  MESH3D  input  data  file  can  be  arbitrarily  divided  into  seven 
groups  as  described  Below.  It  should  be  noted  that  a  — of  750  elements  may  be 
generated  as  a  mesh  without  altering  tha  dimensions  of  tha  arrays  in  the  program. 

Croup  CD  F808LB1  TITLE  i 

e  Inter  problem  title  occupying  not  more  than  80  characters. 

Croup  (2)  C8APS8IAH  QPADlATCkl 

*  Inter  tha  order  of  Cauasian  quadrature  formula  used  for  the  integration  over 
each  element  and  the  corresponding  Gaussian  abscissas  and  coefficients. 

(Rote:  The  Gaussian  abscissae  and  its  coefficients  are  used  in  the  Gaussian 
quadrature  formulae  to  evaluate  various  integrals  over  elnents  and  cells  in 
BIB.  A  short  summery  of  thair  values  is  listed  in  the  Appendix  or  a  compre¬ 
hensive  set  of  tha  coefficients  may  be  found  in  kef  5.  As  a  guide  line  four th- 
or  sixth-order  of  Caaeaian  Quadrature  is  sufficient  for  most  cases.  If  in 
doubt,  the  safest  way  to  get  information  on  the  stability  of  the  results  is  to 
maka  several  runs,  at  the  subsequent  stage  of  using  UX3D,  with  different  order 
of  Gaussian  quadrature.  It  should  be  noted  thet  tha  Gaussian  quadrature  is 
not  used  in  tha  algorithm  of  Che  mash-generation  program  so  they  can  be  neg¬ 
lected  here  by  entering  a  ZOO  instead;  however  it  will  need 'to  be  changed  for 
tha  IH3D  ran.) 

Group  Q)  hatehc.  nartxrr 

*  later  the  Tseng's  modulus  and  Folssou’s  ratio. 

(Hotel  Neither  Toung’s  modulus  nor  Poisson* s  ratio  are  used  in  the  algorithm 
of  tha  mesh-generation  progra;  their  values  will  not  therefore  affect  tha 
mesh  result  But  are  needed  for  ISZ3D.) 

Group  (A)  Wl  BATA 

*  Inter  the  total  number  of  points,  total  maber  of  blocks,  type  of  elements  to 
be  generated  (op  eight}  and  the  dimension  of  the  co-ordinate  system  used 

Co f  two  or  three).  These  data  are  read  in  subroutine  9OL0C  and  the  para¬ 
meters  defined  in  the  program  as< 

■F0H  -  Maber  of  eo-ordiaate  points  defining  the  structural  outline. 

HELM  -  Humber  of  blocks  fr— of  50  blocks  may  be  specified). 
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■BOCK  -  Types  of  el  meant  which  art  permitted  in  eh*  program  arts 

fa)  Quadratic  eight  noda  iaopag— a  trie  elmeents, 

fb)  Llnaar  four  noda  elements, 

(e)  tinner  three-node  triangular  element. 

BDDJE  -  Eumber  of  co-ordinate  diaanaiona  (*/  two  or  three) . 

Croup  f3)  Structural  block  data 

•  .  Intar  tha  block  mmber,  tha  numbers  of  tba  nodal  points  of  tha  block.  (Kota  tha 

numbering  convention  for  a  block  is  aa  followa:  as  tha  surface  is  viewed  along  tha 
inward  animal  from  outside  tba  body,  the  nodal  numbers  are  opacified  in  clockwise 
order,  £s  right-hand-screw  rule.) 

These  par  ana  tars  are  defined  la  tba  program  aat 
WIHEI.  -  Structural  block  number. 

LMDS  -  Structural  block  data.  asCDSCmDOEL.rMDE),  DWDE-I  .RHODE) ,  ts  the  block 
topology  definition. 

Croup  (6)  Co-ordinates  of  the  nodal  points 

•  Enter  tba  point  m»bot  and  tha  three-dimensional  co-ordinate  (z,y,s)  of  tha  modal 
points  of  tba  blocks.  Thaaa  parameters  are  defined  in  tha  program  aa: 

JPOTH  -  point  anther. 

C00ED(JP01H,IDIME)  Cn>S9E*J  .HD1ME)  -  Cz,y,z)  coordinates  of  tbs  points. 

(Kota:  marfim  of  200  points  assy  be  used  to  define  the  structural  outline  without 
altering  the  also  of  tba  array  dimension  in  tbs  progrm.) 

croup  <7)  mgamrnmu 

•  Inter  block  maker,  maker  of  alamanta  in  tha  C  direction  and  maker  of 
eleanuca  in  tha  ij  direction. 

a  Enter  the  corresponding  weighting  factors  for  tbs  C  direction  and  the  weighting 
factors  for  tbs  n  direction. 

Those  are  tba  input  data  necessary  to  subdivide  each  structural  block  into  smaller 
alamanta  and  these  parameters  are  defined  in  the  program  aa: 

EELOC  -  Block  number. 

EDIVX  -  Busbar  of  elements  in  C  direction  into  which  the  block  is  to  be  subdivided. 
BOOT  -  as  EDIVX  but  in  {  direction. 

(Bote:  «*»<«■■  value  of  four  may  be  specified  for  EDIVX  and  EDIVT  without  altering 
the  array  dimension  in  tha  program.) 

MUD-CF-nU 

2.6  Teat  case 

A  flat  plaee  with  a  normal  central  bole  is  used  bare  as  aa  asaple  to  darnnstracs 
tha  operating  procedures  of  tha  mesh  generation  program  MESK3D.  Tha  diaanaiona  of  the 
flat  plats  are  shewn  in  Big  S.  In  order  to  make  tha  problem  as  simple  as  possible  tba 
ijmauj  about  the  r-  and  y- ernes  is  used  so  that  only  a  quarter  of  the  plate  is  modelled 
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in  thi*  tamp  In.  Tb*  initial  Mih  of  tb*  flat  plat*  contain*  a  total  of  10  larg*  »truc- 
tnral  block*  and  a  total  of  32  point*  a*  above  in  Tig  7  and  Tig  6. 


In  order  to  generate  a  a*«b  the  largo  atructural  block*  art  eubdividad  aaing  the 
vaighting  factor*  a*  described  in  Metion  2.3.  k  lifting  of  tb*  input  data  file  for 
tbo  NESB3D  progrnn  i*  liatad  bolow  vith  dateriptlona.  The  data  fila  nan*  i*  daaignatad 
a*  rr.Bkii- 


Pata  fila  ffJU 


D*»crlption« 


THU:  TUX  TUB  WITH  k  CBNTBE  BOLT 


Problen  title 


.23MI9I9 

.661209*0 

.932*6969 

-.23661919 

-.66120940 

-.932469*9 

7000.0  0.3 


32  10  8  3 


.46791393 

.36076138 

.17132449 

.46791393 

.36076158 

.17132449 


1  1  2  3  4  3  6  7  8 

2  5  4  3  9  10  II  12  13 

3  12  II  10  14  13  16  17  18 

4  17  16  13-  19  20  21  22  23 

3  7  6  5  13  12  24  23  26 

6  25  24  12  18  17  23  22  27 

7  7  26  23  28  29  30  I  8 

8  23  27  22  21  20  31  29  28 

9  20  19  13  14  10  32  29  31 

10  1  30  29  32  10  9  3  2 

1  10  0 

2  1  0  .1 

3  I  0  .2 

4  1  .4  .2 
3  1  .8  .2 
6  1  .8  .1 
7  I  .8  0 
6  1  .4  0 

9  .3  0  .2 

10  0  0  .2 

11  .4293  .4293  .2 

12  .8386  .8386  .2 

13  .9233  .8152  .2 

14  0  .3  .2 
13  0  I  .2 

16  .4  1  .2 

17  .8  I  .2 

18  .8132  .9233  .2 

19  0  I  .1 

20  0  I  0 

21  .410 

22  .810 

23  .8  I  .1 

24  .8386  .8386  .1 


....  order  of  Cauaaian 
quadrature  and  coofficiant*. 


....  Toung'a  aodulu* 
and  Toitaon'*  ratio 

....  Total  nunbar  of  point*, 
total  nunbar  of  block*,  alanant 
tppa  «a  8)  and  co-ordinat* 
diaanaiona-. 

....  Block  nabar  and  block 
data  apacifiad  in  tb*  right- 
hand-*crav  rul* 


....  Point  nnbtr  and  tb*  tbraa- 
dinanaional  cartaaian  eo-ordinata 

apatan 


i 

i 

i 
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23  .ISM  .8386  0 

26  .9233  .8132  0 

27  .8132  .9235  0 

28  .6293  .6293  0 

29  0  0  0 

31  0  .5  0 

32  0  0  .1 

1  3  6  ....  Block  (1)  and  weighting 

3.5  3  3.5  factor*  In  the  K  and  n  direction. 

2  3.5  3  1.5 

2  6  1 

1.3  3  3.3  2 
I 


....  Block  (2)  and  weighting 
factor*  in  th*  5  and  n  direction. 


361  ....  etc 

1.5  3  3.3  2 
1 

6  6  3 

1.5  3  3.3  2 
3.3  3  3.5 


5  3  1 

3.5  3  3.3 
1 

6  3  1 

3.5  3  3.5 
1 

7  1  6 
1 

1.5  3  3.5  2 

8  16 
1 

1.5  3  3.3  2 

9  3  1 

3.5  3  3.5 

10  1  3 
I 

3.5  3  3.5 

m>  ....  *nd-of-file 


Thi*  data  file  (TP  .DAI)  aay  then  be  aubaitted  to  tha  MESH3D  prograa  u*ing  eh*  eomaod*  a* 
follow*! 


laaaign  TP .DAI  for 003 
8aa*ign  TP. MBS  for 006 
Iron  MESH3D 


(a**ign  to  logical  unit  3) 
(aaiign  to  logical  unit  6) 
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The  input-data  C77.D4T)  will  ba  checked  for  tha  total  umber  of  nodal  point*  and 
eleaanca  to  ensure  tba  array  dimensional  aixa*  in  tha  MESH3S  prograa  will  not  ba 
axcaadad.  da  excessive  vain*  of  tba*a  umbers  will  causa  an  execution  arror  and  an 
arror  aassage  will  Ba  displayed  before  aborting  tba  prograa.  If  tba  total  noaher  of 
nodal  points  and  alaaant*  ia  not  too  large  then  tba  aaab  coordinate*  of  tba  three- 
diaanaional  aaab  will  Ba  autoaatically  generated  according  to  tba  process  of  sub¬ 
division  and  tba  corresponding  weighting  factors.  The  output  data  file  containing  tha 
generated  wash  data  of  tha  quarter  flat  plate  ia  called  (FP.MES) .  Before  this  file  can 
be  snBnieted  for  the  BXK3D  run  there  are  two  aore  steps  the  user  any  need  to  carry  out: 

0)  visual  inspection  of  tba  topology  of  tba  generated  wash,  and 
(2)  specify  Boundary  conditions. 

Tba  visual  inspection  of  tba  aeah  can  ba  dona  once  tha  aasb  is  plotted  on  paper.  The 
procedures  for  cha  plotting  of  tba  aaab  on  paper  are  described  in  section  3.  Once  the 
aasb  baa  Been  checked  for  aiaaligunant,  anoaalous  points  or  any  incaapatibility  of 
alaaant  distribution*  ate  it  can  than  be  aodified  to  include  tba  boundary  conditions  of 
tba  profiles;  tba  Boundary  conditions  are  inserted  according  to  the  f onset  specified  for 
boundary  conditions  in  tha  BXE3S  users'  guide1 .  Finally  the  user  can  suboit  tba  aodified 
QrP.MES)  input  data  file,  OfflDTP.HES)  say,  and  run  tba  B1Z3S  prograa  using  the  comends 
as  follows: 

^assigns  HODFF.MBS  forOOS  (assign  to  logical  unit  5) 

lasaign  TP.BES  for006  (assign  to  logical  unit  6) 

ntOK  BTE3D 

where  ’FF.BES'  is  the  output-data-flle  of  tba  BIE3S  prograa  containing  the  displscments, 
tha  tractions  and  the  surface  stresses  for  all  tha  nodes  of  tha  threa-diaensional  aesh 
of  the  quarter-plate. 

3  MESH  PLOTTING 

The  checking  of  tha  nodal  co-ordinates  of  three— disantioual  aesh  data  1*  often 
tedious  and  difficult.  Ia  particular,  if  tba  aesh  data  ba*  been  generated  by  a  aesh 
generation  prograa,  the  large  voluae  and  tba  three-diaausional  coupler! ty  of  the  aesh 
data  would  aaka  tha  checking  process  an  awescae  task.  The  obvious  alternative  approach 
i*  to  plot  the  generated  aesh  on  paper  ia  an  isoastric  view.  Tbs  visual  topology  of 
the  aaab  would  allow  for  tha  checking  of  aisaligmsnt,  anoaalous  point*  or  any  incom¬ 
patibility  of  alaaant  distributions  ate.  Such  a  visual  presentation  for  a  aasb  in 
three  diaension*  is  invaluable  and  absolutely  essential  for  any  aesh  design  and  aasb 
generation  work.  Indeed,  during  tha  plotting  process,  if  one  follows  tba  path  of  the 
plotting  pen  on  the  chart,  on*  can  easily  datsraiaa  without  recourse  to  tha  nmarical 
data  file,  whethar  eha  aesh  ha*  Been  generated  ia  the  specified  ouabering  convention. 
Vurtfiaxaorc,  with  tha  option  of  having  all  the  nodal  nunbert  of  all  the  alaaant  ambers 
or  both  plotted,  the  nmberad  aasb  can  ba  used  to  assist  tbs  specification  of  all  the 
boundary  conditions  in  which  each  individual  nodal  point  or  surface  slaaent  can  be 
precisely  located  and  constrained  if  required. 
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A  three-dimensional  neah  plotting  coaputar  prograa  (DKAH3D)  bu  ban  written  and 
implemented  is  tba  local  coaputar  system.  Tba  prograa  ia  written  la  BASIC  and  It  rant 
interactively  on  a  Hewlett-Packard  D2647F  graphics  terminal.  It  can  either  ha  need  to 
plot  a  tbree-diaanaional  mash,  in  an  laoaatric  viaw,  on  a  graphic  terminal  diaplay  or 
on  paper,  in  thia  a action,  tba  handling  of  input  and  output  data  filoa  mill  be  explained 
and  horn  exactly  the  graphic  plotting  pro grama  can  be  need  in  conjunction  mi th  the  aaah 
gonaracion  prograa  QS3R3S)  and  the  boundary  integral  aquation. prograa  (B1X3D). 

3*1  The  graphic  plottlnx  coaoutar  prograa  (DBAM3D) 

The  coaputar  prograa  BBAB3D  can  be  uaad,  on  a  Hemlett  Packard  BP2647F  terminal,  to 
plot  a  three-dimensional  viaw  of  a  aaah.  The  plotting  prograa  0DHAH3D)  ia  nom  available 
in  the  local  coaputar  ayeten.  It  haa  been  axtenaivaly  uaad  throughout  the  preaent  work 
and  ahomn  to  be  maaful  tor  checking  of  a  generated  aeah.  The  prograa  can  be  uaed  inter¬ 
actively  By  the  uaera  ao  that  tha  view  anglaa  of  an  iaoaaerie  plot  can  be  rotated 
through.  0  to  360  Segreea  in  any  axe a  to  a  horn  any  aide  of  tha  figure.  41  ao,  the  diman- 
aiona  of  tha  figure  in  tha  X-  and  y-planaa  can  ba  diatorted  ao  that  exaggeration  of  any 
particular  plana  of  the  figure  nay  be  aede.  Phan  tha  final  projection  of  the  configura¬ 
tion  ia  accepted,  it  can  Ba  plotted  on  paper  aa  a  permanent  record. 

3.2  Bata  preparation  for  the  plotting  prograa 

A  BIX3D  input  date  file  r aquiree  very  little  change  before  it  can  be  aubaittad 
for  a  plotting  prograa  run.  In  fact  the  only  change  which  ia  required  ia  to  convert 
a  TCKTBAM  input  data  file  fornat  into  a  BASIC  input  date  file  foraat.  The  converaion 
ia  an  aaaential  atap  and  it  ia  dona  By  aeparating  each  data  item  by  a  coana. 

The  following  section  mill  liat  all  the  atepa  required  to  prepare  a  plotting  prograa 
input  data  file. 

3.3  Operating  procedurea  for  the  graphic  plottlnx  program 

(1)  Separate  the  aaah  data  in  a  B1X3D  input  data  file  by  c cease .  Thia  can  be 
achieved  either  by  typing  the  ccaaaa  in  by  hand  or  aubaitting  tha  input  data  file  to 
a  utilities  prograa  called  BBADCOB  ia  tba  TAX  coaputar.  The  utilities  progra  will 
separata  all  the  data  in  the  data  file  with  coenss  automatically: 

las  sign  (BZB3D  input  data  file  )  WB001 

♦run  BXABCOB 

and  the  Output  of  the  aodified  aeah  data  ia  assigned  to  the  display. 

(2)  Mom  me  can  transfer  thia  aodified  aaah  data  frea  tha  diaplay  to  a  floppy  disc, 

(3)  and  then  run  the  BASIC  plotting  prograa  by  typing'  in 

>  GtT  "BBAB3D"  (asauaing  that  tha  plotting  progra  ia  already  ia  the  disc) 

>  BOX 

At  this  point,  the  plotting  prograa  will  respond  with  prompts  la  an  intaractlve  node 
and  a  aeries  of  prompts  will  be  displayed  on  the  screen: 
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T 
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>  Plotting  device?  (diep lay-0  or  plotter-5) 

>  Input  file  umi 

>  Slant  type?  (Triangular-3  or  Quadrilateral-8) 

One*  thooo  quoationo  ora  aaaaarad,  the  progran  will  taka  over  and  plot  out  tha  three- 
dinaneional  aiaw  of  pour  naah.  After  It  haa  finished,  further  profit a  will  appaar  on 
tha  display: 


>  Do  you  wish  to  have  tha  elment-nunber  laballadt  Ct  or  8) 
*  Do  you  viah  to  have  tha  nodal -ounber  laballadt  (T  or  8) 


Oban  thaaa  ara  dona  tha  nesh  plotting  proeaaa  ia  eoapletad  and  tha  progran  will  r a-* cart 
with  tha  following  pronpta: 

>  DIES  view-angles  (rotation  about  a*  and  y-axaa)  and  diatortion  ia  y  diraction: 

>  ESTER  Plotting  device?  (displey-0  or  plottar-5) 

>  Input  total  number  of  hiddan  alaaianta: 

>  Do  you  want  to  plot  diatortad  ahapaf  (T  or  8) 

The  uaar  can  aalact  any  ona  of  thaaa  optiona  and  carry  on  aa  bafora. 


3.4  Taat  caaa 

Tha  fila  callad  (FP.MES)  ia  tha  ganaratad  aaah  data  file  for  the  quarter-plate 
azaapla  aa  daacribed  in  aaction  2.6.  It  ia  uaad  hare  to  danonatrsta  the  operating  pro- 
cadurea  of  tha  plotting  computer  progran  DIA83D. 

Bafora  tha  data  fila  can  be  aubnitted  for  tha  plotting  run,  the  fomat  of  the  data 
fila  nuat  be  converted  to  the  BASIC  data  fomat  aa  daacribed  in  aaction  3.3  The  converted 
data  of  fila  (TP.ME3)  ia  liated  below:  it  ahould  be  noted  that  thia  Hating  baa  baan 
coudenaad  to  aava  apace. 


TOo  threo-diaenaional  praaantationa  of  tha  ganaratad  naah  are  plotted  on  paper  aa 
abowu  in  Pig  9  and  Pig  10.  .It  can  be  eeen  that  a  total  of  52  quadrilateral  element «  and 
a  total  of  1  SB  nodal  pointa  have  been  ganaratad  ia  thia  caaa.  The  element  mabering 
aye ten  end  nodal  naaberiag  ayatan  of  tha  generated  naah  ara  eboon  ia  Piga  9  and  10, 
raapactivaly. 


Data  fila 


Daecrlptione 


TITLE:  PUT  PLATE  WITH  A  CESTRE  BOLE 


....  Title 


138.  52 


....  Busbar  of  nodaa  and  mmber  of 
elananta. 


1,  1.000000 

2,  1.000000 

3,  1 .000000 


0.0000000*00, 

0.00000001*00, 

0.00000001*00, 


0.00000008*00 
3. 50 00 00 0E -02 
7.00000001-02 


....  Bode  unbar  and 
a.y.a  cartesian  co-ordinate a 


155,  0.000000 

156,  0.000000 


0.00000OOE*O0,  0.1300000 
0.00000008*00,  0.1630000 


u 


157. 

0.500000 

§ 

0.00000008*00  ,  7.00000008-02 

158. 

0.500000 

• 

0.00000008*00  ,  0.1300000 

1. 

1. 

2. 

3. 

9,  14.  13,  12,  8 

3. 

3, 

*. 

-  3. 

10,  16,  IS,  34,  9 

3. 

5. 

6. 

7. 

11,  18,  17,  16,  10 

4. 

12. 

13. 

14. 

20,  25  ,  24,  23,  19 

....  Blaaant  sober  od  alwcnc 
data  specified  la  right-hand- 
scree  conrantion. 


48,  107  ,  98.  93,  1S1,  155.  154.  353.  150 


*9. 

93,  84,  79,  70,  65,  136,  135, 

15J 

30. 

1,  IU,  145,  152,  153,'  157, 

3. 

2 

31. 

3,  157,  153,  154,  155,  158, 

3. 

4 

32. 

33. 

8BD 

5,  158,  155,  156,  63,  36, 
20,  1 

7, 

6 

«...  ritm  agios  ad  distortion. 

•pea  fad  Of  fill. 

4  cogarogg  mm 

Tha  instructions  contained  la  this  HaanrandiM  will  aaabla  uaars  to  xaa  the  progress 
MESH3D  and  0UM3D  aa  they  ataoB,  to  ohtala  respectively  a  aeah  aad  a  graphic ■  prssenta- 
tioa  of  a  aeah  of  a  staple  or  coaplex  threa-diaeaaloaal  configuration. 

A  throa-dlaaaaloaal  aash  pros so  tat Ion  of  a  real  structure  eaa  ba  vary  eoaplan 
and  difficult  to  haadla  aaanally,  hat  by  uaing  those  prograas  tha  usars  can  leave  tha 
nahar  aanipnlatlon  to  tha  prograas  aad  concentrate  on  tha  accuracy  of  tha  asdalliag  of 
tha  problsa. 


1 


1H  HS  IU/ / 


Ciuilu 

abtcltua 


Staulu  coefficient* 


4 

0.3399810435 

0.8811363115 

-0.3399810435 

-0.86113631)5 

0.6521451548 

0.3478548451 

0.6521451548 

0.3478548451 

6 

0.2386191860 

0.6612093864 

0.9324695142 

-0.2386191860 

-0.6612093864 

-0.9324695142 

0.4679139345 

0.3607615730 

0.1713244923 

0.4679139345 

0.3607615730 

0.1713244923 

8 

0.1834346424 

0.5255324099 

0.7966664774 

0.9602898564 

-0.1834346424 

-0.5255324099 

-0.7966664774 

-0.9602898564 

0.3626837833 

0.3137066458 

0.2223810344 

0.1012285362 

0.3626837833 

0.3137066458 

0.2223810344 

0.1012285362 
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Fig  8  Preliminary  mesh  design  and  nodal  numbering  system  of 
a  quarter  of  the  flat  plate 
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fig  10 


FI9  10  Mesh  design  and  nodal  numbering  system  of  a  ouarter  of 
the  flat  plate 
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